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Abstract 

We explore the thermodynamics of dark energy taking into account the existence of the observer's 
event horizon in accelerated universes. Except for the initial stage of Chaplygin gas dominated 
expansion, the generalized second law of gravitational thermodynamics is fulfilled and the temper- 
ature of the phantom fluid results positive. This substantially extends the work of Pollock and 
Singh on the thermodynamics of super-inflationary expansion. 
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I. INTRODUCTION 



Nowadays there is an ample consensus on the observational side that the Universe is un- 
dergoing an accelerated expansion likely driven by some unknown fluid (called dark energy) 
with the distinguishing feature of violating the strong energy condition, p + 3P > [lj. The 
strength of this acceleration is presently a matter of debate mainly because it depends on 
the theoretical model employed when interpreting the data. While most model independent 
analyses suggest it to be below the de Sitter value [J| it is certainly true that the body of 
observational data allows for a wide parameter space compatible with an acceleration larger 
than de Sitter's 3,0]. If eventually this proves to be the case, the fluid driving the expan- 
sion would violate not only the strong energy condition but the dominant energy condition, 
p + P > 0, as well. In general, fluids of such characteristics, dubbed "phantom" fluids, face 
some theoretical difficulties on their own such as quantum instabilities 15]. Nevertheless, 
phantom models have attracted much interest, among other things because in their barest 
formulation they predict the Universe to end from an infinite expansion in a finite time ( "big 
rip"), rather than in a big crunch, preceded by the ripping apart of all bound systems, from 

. While this scenario might look weird, given our 
incomplete understanding of the physics below p + P = and the scarcity of reliable obser- 
vational data, it should not be discarded right away; on the contrary we believe it warrants 
some further consideration. 

Recently, it has been demonstrated that if the expansion of the Universe is dominated by 
phantom fluid, black holes will decrease their mass and eventually disappear altogether [7[. 
At first sight this means a threat for the second law of thermodynamics as these collapsed 
objects are the most entropic entities of our world. This short consideration spurs us to 
explore the thermodynamic consequences of phantom-dominated Universes. In doing so 
one must take into account that ever accelerating universes have a future event horizon (or 
cosmological horizon). Since the horizon implies a classically unsurmountable barrier to our 
ability to see what lies beyond it, it is natural to attach an entropy to the horizon size (i.e., 
to its area) for it is a measure of our ignorance about what is going on in the other side. 
However, this has been proved in a rigorous manner for the Sitter horizon only [8] which, in 
addition, has a temperature proportional to the Hubble expansion rate, Tu = fiHlil-Rkn). 

nnnn 

Nevertheless, following previous authors |y, HO, IUJ here we conjecture that this is also 
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true for non-stationary event horizons. While certainly this is a bold assumption that we 
cannot justify at present we believe the reasonableness of the results lend support to it. 
Phantom expansions of pole-like type 



a(t) oc — (t < £*, < n = constant) , (1) 

as proposed by Pollock and Singh and Caldwell jl^, arise when the equation of state 
parameter w = P 7 p = constant < —1. Current cosmo logical observations hint that w 

n 

may be as lower as —1.5 |13J. In what follows it will be useful to bear in mind that n = 
-2/[3(l+w)\ > 0. 

From the above equation for the scale factor it is seen that the Hubble expansion rate 
augments 

whereby the radius of the observer's event horizon 

Re = <t) = a(t) P*j- = ±ll = -^-H^ (c = 1) (3) 

Jt a(t') J t a(t') 1 + n l + n 

decreases with time, i.e., Rh < 0, and vanishes altogether at the big rip time t*. 
Consequently the horizon entropy, 

». = 7F' (4) 

where A = AttR 2 h denotes the area of the horizon and £ p i = a/ hG / c 3 is the Planck's 
length, diminishes with time, Sh < 0. This is only natural since for spatially flat FRW 
phantom-dominated universes one has 

H = -^k( P + P)>0 (k = 8tvG), 

as these fluids violate the dominant energy condition. Then the question arises, "will the 
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generalized second law (GSL) of gravitational thermodynamics, S + Sh > 0, be satisfied?" 
The target of this work is to answer this question. According to the GSL the entropy of 
matter and fields inside the horizon plus the entropy of the event horizon cannot decrease 
with time. To our knowledge the first authors who considered this question assu ming energy 
sources that violate the dominant energy condition were Pollock and Singh ldl. They 
studied the super-inflation models of Starobinsky [3| and Shafi and Wetterich jl5|| in the 
regime that they departed only slightly from de Sitter expansion and found that while 
the former fulfills the GSL the latter does not. We extend the work of these authors by 
considering other models of accelerated expansion (phantom and non-pantom) not restricting 
ourselves to small deviations from de Sitter. As it turns out, except for the Chaplygin gas 
dominated universe (section III) the GSL is fulfilled in all the instances explored. 

Since sooner or later the expansion will get dominated by the dark energy fluid we will 
neglect, for the sake of transparency and simplicity, all other sources of energy (e.g., matter, 



radiation, and so on). For the interplay 

nnnn 

the cosmological event horizon, see Refs. ^laimiii 



je ween ordinary matter (and/or radiation) and 



II. DARK ENERGY WITH CONSTANT w 

We begin by considering a phantom fluid inside the cosmological event horizon of a 
comoving observer. Its entropy can be related to its energy and pressure in the horizon by 
Gibbs' equation 

TdS = dE + Pd(lnR 3 H ) . (5) 

In this expression, owing to the fact that the number of "phantom particles" inside the 
horizon is not conserved, we have set the chemical potential to zero. From the relation 
E = (4n/3)pRj { , together with the Friedmann equation p = (3/k)H 2 , Eq. (jHJ), and the 
equation of state P = wp (with w < —1), we get 

TdS = -—^—dR H . (6) 
k 1 + n 

Since dR H < the phantom entropy increases with expansion (so long as T > 0). 
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One may argue that it is wrong to apply Gibbs equation to a fluid that in reality is a 
phenomenological representation of a scalar field assumed to be in a specific state, and that 
in such a case its entropy should strictly vanish. This reasoning could be correct if we were 
supposing that the fluid phenomeno logically represents just one single field and this one in 
a pure state, but this is far from being the more natural assumption. We see the fluid as a 
phenomenological representation of a mixture of fields, each of which may or may not be in 
a pure state but the overall (or effective) "field" is certainly in a mixed state and therefore 



entitled to an entropy. This is the case, for instance, of assisted inflation [17 1. 

To proceed further, we must specify the temperature of the phantom fluid. The only 
temperature scale we have at our disposal is the temperature of the event horizon, which we 
assume to be given by its de Sitter expression [8] 

T = J± * L (7) 

H 2vr k B 1 + n R H ' 

though in our case H > 0. Thus, it is natural to suppose that T cx Th and then figure out 
the proportionality constant. As we shall see below, this choice is backed by the realization 
that it is in keeping with the holographic principle Q|. As in Ref.Q], the simplest choice 
is to take the proportionality constant as unity which means thermal equilibrium with the 
event horizon, T = Th- In general, two systems must interact for some length of time 
before they can attain thermal equilibrium. In the case at hand, the interaction certainly 
exists as any variation in the energy density and/or pressure of the fluid will automatically 
induce a modification of the horizon radius via Einstein's equations. Moreover if T ^ Th-, 
then energy would spontaneously flow between the horizon and the fluid (or viceversa), 
something at variance with the FRW geometry. 

After integrating Eq. (jUJ) and bearing in mind that S — > as t — > £*, the entropy of the 
phantom fluid can be written as S = —k B A/ (4 

Some consequences follow: (i 
already noted by other authors 



The phantom entropy is a negative quantity, something 
10, Il9j . and equals to minus the entropy of a black hole 
of radius Rh- (ii) It bears no explicit dependence on w and, but for the sign, it exactly 
coincides with the entropy of the cosmological event horizon. (Hi) Since —S cx E 2 it is not 
an extensive quantity. Note that S does not scale with the volume of the horizon but with 
its area. One may argue that the first consequence was to be expected since (for T > 0) it 



readily follows from Euler's equation T s = p + P, s being the entropy density. However, 
it is very doubtful that Euler's equation holds for phantom fluids since it is based on the 
extensive character of the entropy of the system under consideration (^(J and we have just 
argued that for phantom fluids this is not the case. 

Two further consequences are as follows: (iv) The sum S + Sh vanishes at any time, 
therefore the GSL is not violated -the increase of the (negative) phantom entropy exactly 
offsets the entropy decline of the event horizon, (v) \ S | saturates the bound imposed by 
the holographic principle The latter asserts that the total entropy of a system enclosed 
by a spherical surface cannot be larger than a quarter of the area of that surface measured 
in Planck units. (For papers dealing with the holographic principle in relation with dark 



energy, see 



2l\ and references therein). In this connection, it is interesting to see that if the 



equation of state were such that the entropy obeyed S* = xS, with x a positive-definite 
constant, then the GSL would impose 1 < x. This leads us to conjecture that the entropy 
of phantom energy is not bounded from above but from below, being — Sh its lower limit. 

Nevertheless, in a less idealized cosmology one should consider the presence of other forms 
of energy, in particular of black holes and the decrease in entropy of these objects |7]. It 
is unclear whether in such scenario the GSL would still hold its ground. Nonetheless, one 
may take the view that the GSL may impose an upper bound to the entropy stored in the 
form of black holes. At any rate, the calculation would be much more involved and lies 
beyond the scope of this paper. Among other things, one should take into account that the 
scale factor would not obey such a simple expansion law as and that the black holes 
would be evaporating via Hawking radiation which would also modify the expansion rate 
and, accordingly, the horizon size. 

A straightforward and parallel study for a non-phantom dark energy-dominated universe 
with constant parameter of state (lying in the range —1 < w < —1/3), shows that in this 
case S + Sh = as well [34]. There are two main differences, however; on the one hand the 
entropy of the fluid decreases while the area of event horizon augments, and on the other 
hand S > 0. The latter comes from adopting the view that the dark energy must vanish 
for T — > (Planck's statement of the third law of thermodynamics) and realizing that this 
happens for t — > oo. Again, S saturates the bound imposed by the holographic principle. 

Obviously, one might adopt the view that the phantom temperature should be negative 
But, as mentioned above, this would destroy the FRW geometry. On the other hand, 
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negative temperatures are linked to condensed matter systems whose energy spectrum is 
bounded from above and may therefore exhibit the phenomenon of population inversion, 
i.e., their upper energy states can be found more populated than their lower energy states 



24j . However, all models of phantom energy proposed so far assume some or other type of 



scalar field with no upper bound on their energy spectrum. In addition, while population 
inversion is a rather transient phenomenon the phantom regime is supposed to last for 
cosmo logical times. Moreover, bearing in mind that w < — 1 and dRn < 0, it follows from 
Eq. (jUJ) that if T were negative, then the phantom entropy would decrease with expansion 
and the GSL would be violated. 



III. DARK ENERGY WITH VARIABLE w 

One may argue that the result S + Sh = critically depends on the particular choice 
of the equation of state parameter w to the point that if it were not constant, then the 
GSL could be violated. Here we shall explore this issue taking the phantom model with w 
dependent on time of Sami and Toporensky 2^1 . 

In this model a scalar field, <fr(t), with negative kinetic energy, minimally coupled to 
gravity sources with energy density and pressure given by 

, = -^ + v<#), p = -^-v W , 

respectively, is adopted. The equation of state parameter 

_ P _ V{<t>) 

is lower than —1 as long as p is positive and depends on time. 



When the kinetic energy term is subdominant ("slow climb" 25] ) the evolution equation 
of the phantom field simplifies to 

3Ht(t) - (8) 



Assuming also the power law potential V(4>) = V <p a with a a constant, restricted to 
the interval < a < 4 to evade the big rip, the equation of state parameter reduces to 
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w ~ — 1 — sj, where x = Kip 2 /{2a) is a dimensionless variable, and the field is an ever 
increasing function of time, namely, 



0(t) = 
<p(t) = ^ exp 



y + jg ( 4 -" )a ( t - ti) 



3k 2 
V o/t 



(0 < a < 4), 



(a = 4), 



(9) 
(10) 



where the subscript i stands for the time at which the phantom energy begins to dominate 
the expansion. Thus w increases with expansion up to the asymptotic value —1. Likewise 
the scale factor obeys 



a(4>) = a,i exp 



2a 



As it should, the Hubble factor 



(11) 



H 



2Vq /K\|-f 



3 V2 



a a/4 x a/4 



(12) 



augments with expansion while the cosmological event horizon 



(13) 



decreases monotonically to vanish asymptotically -see Figure □ Here T is the 

incomplete Gamma function 



As in the preceding section, the entropy of this phantom fluid can be obtained by using 
Gibbs's equation and T = Th = (ch/2nkB)H. Again, we get a negative quantity that 
increases with time, namely, 



s^JJ^c* fV(l 

hn J x \ 



. dC I 6 dec . 

4' J 2x 



(14) 



here C = ^3/(21^) (fi:/2)-| + t a~ a 



n 
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FIG. 1: Radius of the event horizon vs x for a = 1, 2, 3 and 4, from top to bottom. In all the 
cases Rh — > as t — ► oo. For convenience sake we have set Xi = 1. 

From the above expressions, it follows that 



r(i^,,)e-(2 + ^)-2^ 



i2 H i- (15) 



Since the quantity in square brackets is positive-definite for any finite x and x > 0, the GSL 
in the form S + Sh > is satisfied. The equality sign occurs just for t — » oo, i.e., when R H 
vanishes. 

As figure 121 illustrates we have that | S |> Sh] however, strictly speaking the holographic 
principle is respected as 5 is always lower than Sh- 

As an example of non-phantom dark energy with variable w we consider the Chaplygin 
gas 2t| . This fluid is characterized by the equation of state P = —A/ p. Integration of the 
energy conservation equation yields p = ■sj A + (B/a 6 ), whereby w = —A a 6 /(A a 6 + B). 
Here a is the scale factor normalized to its present value and A and B denote positive-definite 
constants. Accordingly, this fluid interpolates between cold matter, for large redshift, and a 
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FIG. 2: Evolution of the ratio between the absolute value of the phantom entropy and the horizon 



entropy for the model of Ref. |2|J. The graphs correspond successively to a values 1, 2, 3, and 
4 from bottom to top. For convenience sake we have set Xi = 1 in drawing the graphs. All the 
curves asymptote to x = 1. 



cosmological constant, for small redshift. This has led to propose it as a candidate to unify 



dark matter and dark energy 



28] 



The radius of the cosmological event horizon reads 



R 



H 



3 A 1 / 4 



x 



1/6 



1.4712 



x 



1/12 



1 1 12' 4' 12' 



-x 



(16) 



where 2 -Pi is the hypergeometric function and x = Aa 6 /B is a dimensionless variable. As 
figure H3 illustrates, Lim^oo Rh = \/3/k A -1 / 4 which is nothing but the value taken by H~ l 
in that limit. This was to be expected: for large scale factor the Chaplygin gas expansion 
goes over de Sitter's. 



The horizon temperature is 



10 
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0.9BS 



10 20 30 40 50 60 70 



FIG. 3: Evolution of the adimensionalized radius of the event horizon in the Chaplygin gas model. 
It tends to unity for x — > oo. 



H 



2ir k B V 3 1 x 



(17) 



while the entropy of the fluid and the entropy of the horizon evolve as 



dS 



A8k B ix 2 1 
'h^A 1 / 2 x 2 /*(l + x) 



1.4712 2 F, [-,-,-; -x 



dx, (18) 



and 



dS 



H 



l6k B ir 2 
hn 2 A 1 / 2 



X 



1/6 



1.4712 -x 1/2 2 Fi 



x 



2x" 5 / 6 



/ll 13 

1.4712-^ 2Fl (___._ x 



1 1 13 
12'4' 12' _X 



a;3/4( 1 + a .)l/4 j ' 



(19) 



respectively. The two latter equation imply that S + S# > for x > 2.509 or equivalently, 
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for pi < 1.1 8 A 1 / 2 . Thus, there is an early period in which the GSL seems to fail. However, 
this might be seen not as a failure of the GSL but as an upper bound on the initial value of 
the energy density of the Chaplygin gas. 

Notice that, S is a positive, ever-decreasing function of the scale factor that complies 
with the holographic bound S < Sh- 



IV. THE ISSUE OF NEGATIVE ENTROPIES 



The fact that S < looks troublesome. In particular, it tells us that the entropy of 
the phantom fluid is not to be understood as a measure of the number of microstates 
associated to the macroscopic state of the fluid, i.e., the well-known statistical mechanics 
formula S = fee In ll 7 breaks down. However, the fact remains that the entropy of systems 



violating that condition (w < —1) is found to be negative (see also, [lCj, |l9j|) whence one 
is rightly entitled to wonder whether these systems can be realized in nature. While the 
answer to this question lies beyond our present capabilities we may try to advance some 
ideas, pending a deeper study. First, we wish to recall that no physical law asserts that 
the entropy of a system has to be positive, the second law of thermodynamics simply 
establishes that the entropy of isolated systems cannot decrease. Secondly, one may draw 
a rough formal analogy between dark energy fluids with equation of state P = wp and a 
very well known physical system, the mono-atomic ideal gas. The entropy of the latter 
is given by the Sackur-Tetrode equation (see Eq. (9.54) of Ref. which can be written as 



k B 2 



2nK 2 \ 3/2 

n 



where a is the entropy per particle, m the particle mass, and n the particle number density. 
Clearly, by lowering the temperature below T c = (2 7rh 2 /em) n 2 ^ 3 the entropy becomes 
negative and the pressure, P = 2ne/3, with e the energy per particle, decreases accordingly 
(though it stays positive). Thus, by lowering P (at fixed n) negative entropies can be 
formally attained. Obviously, one can argue that before a negative entropy state could be 
reached quantum effects would invalidate the approximations leading to the Sackur-Tetrode 
equation, and that, in any case, the gas would condensate ahead of that state. 

This is somewhat analogous to a dark energy fluid whose pressure is steadily decreased 
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by lowering w. When the w = —1 divide is crossed, the entropy becomes negative but 
simultaneously quantum instabilities might destroy the state 

So, we are at a cross-road: either states violating the dominant energy condition are to 
be banned (as they entail negative entropies) or statistical mechanics is to be generalized to 
encompass these states. It is for the reader to decide which path to follow. 

V. QUASI-DUALITY BETWEEN PHANTOM AND NON-PHANTOM THER- 
MODYNAMICS 

The duality transformation 



of spatial FRW universes invariant and the scale factor transforms as a ^ a = 1/a [30]. 
Thus, a phantom dominated universe with scale factor a = l/(—t) n (for simplicity we 
assume a phantom scale factor as Eq. (0) with t* — and t < that tends to 0~) becomes 
a contracting universe whose scale factor obeys a = {—t) n . This universe begins contracting 
from an infinite expansion at the infinite past and collapses to a vanishing scale factor for 
t = 0~ (left panel of Fig. gj. 

Reversing the direction of time (i.e., the future of t, t + dt, transforms into t — dt), 
the contracting universe is mapped into an expanding one whose scale factor grows with 
no bound as t tends to — oo. Notice that the first transformation takes the universe from 
expanding to contracting and the second transformation makes the universe expand again, 
i.e., the final cosmology has H > 0. It should be remarked that the latter transformation 
also preserves Eqs. (J2H)- 

Under these two successive operations, the final and the initial equation of states pa- 
rameters are related by w — > W = — (w + 2). Consequently, phantom dominated universes 
with n > 1, (i.e., —5/3 < w < — 1) are mapped into non-phantom dark energy dominated 




(20) 



leaves the Einstein's equations 



H 2 = IP, 

p + 3H(p + P) = 0, 



(21) 
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universes (see Fig. 0}. Therefore, both of them, the original and the transformed universes, 
have an event horizon and entropies whose relation we derive below. 

The radius of the event horizon of the original universe is Rh = {—t)/ (1 + n) meanwhile 
the radius of the event horizon of the transformed universe is Rh = (—t) /(n — 1) = Rh[(1 + 
n)/(l — n)]. The transformations preserve the horizon temperature as the Hubble factor 
does not change. From the definition of the horizon, it is readily seen that 

^ = (22) 

In virtue of the above expression together with the Gibbs equation one obtains 

T f = " 47r(P + P) ^' (23) 
whereby, the entropy transforms as 

S^S= -[(l + 3w)/{5 + 3w)] 2 S . (24) 



Therefore, the entropy of the final non-phantom dark energy dominated universe is propor- 
tional, but of the opposite the sign, to the entropy of the original phantom universe. We 
might say that the duality transformation "quasi" preserves the thermodynamics of dark 
energy. This result is in keeping with the findings of Section |H] 



VI. CONCLUDING REMARKS 

We have substantially extended the study of Pollock and Singh |l(Lwho considered 
the models of super-inflation (characterized by H > 0) of Starobinsky [14j and Shafi and 
Wetterich Q], by taking up models of accelerated expansion (both with H > and with 
H < 0) with w constant as well as models with w variable (such as the model of Sami and 

n n 

Toporensky and the Chaplygin gas model [27]). 

We have found that, irrespective of whether w is constant or not, phantom fluids (as- 
sociated with H > 0) possess negative entropy, transcend the holographic bound in the 
sense that | S |> Sh and their temperature and entropy increase as the Universe expands. 
By contrast, non-phantom fluids (associated with H < 0) have positive entropy, satisfy 
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t t+dt t t+dt 
> < 

FIG. 4: The left panel depicts the scale factor of the original phantom dominated universe (solid 
line). The dashed line shows the evolution of the scale factor dual of this phantom DE model, a 
contracting universe. In both models (the phantom and its dual) the time grows from left to right 
as tends to t = 0~ and the instant t + dt lies at the future of the instant t. In the right panel, 
the time direction is reversed (t + dt lies at the past of t). Under this operation, the contracting 
model of the left panel becomes the non-phantom DE dominated universe of the right panel if 
-5/3 < w < -1. 

the holographic bound, S < Sh, and their temperature and entropy decrease with expan- 
sion. Except for the Chaplygin gas model which violates the GSL at the earliest stage of 
dominating the expansion, the GSL is fulfilled in all the cases. 

It goes without saying that negative entropies are hard to assimilate, they have no clear 
physical meaning in the context of statistical mechanics; especially because the Einstein- 
Boltzmann interpretation of entropy as a measure of the probability breaks down. Systems 
of negative entropy appear to lie outside the province of statistical mechanics as is currently 



formulated [2^, |29|, |3l|. However, if future cosmic observations conclude that definitely 
w < —1, it will become mandatory to generalize that subject accordingly. We wish to 
emphasize that the laws of thermodynamics do not entail by themselves that S ought to be 
a positive quantity. The latter will be found to be positive or negative only after combining 
these laws with the equation (or equations) of state of the system under consideration. What 
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we believe to be at the core of thermodynamics is the law that forbids the entropy of isolated 
systems to decrease, this law encompasses the GSL for gravitating systems with a horizon. 
Another interesting result is that, unlike recent claims [23 , the temperature of the phantom 
fluid must be positive if the GSL is to be satisfied (see Eq. (jUJ)). In our view, this settles 
the dispute between those rejecting phantom energy on grounds that it cannot be physically 



acceptable because its entropy is negative 19 and those who evade negative entropies by 
advocating negative temperatures instead 23j]. Both approaches overlook the role of the 
event horizon. 

As pointed out by Nojiri and Odintsov, it may well happen that in the last stages of 
phantom dominated universes the scalar curvature grows enough for quantum effects to 
play a non-negligible role before the big rip. If so, the latter may be evaded or at least 
softened j^. While our study does not incorporate such effects they should not essentially 
alter our conclusions so long as the cosmic horizon persists. 

It is noteworthy that the duality transformation p + P — > — (p + P) -see, e.g. j^oj]- along 
with reversing the direction of time (i.e., dt — ► —dt) leaves Einstein's equations invariant and 
maps phantom cosmologies, with both H and H positive, into non-phantom cosmologies with 
H > but H < . However, this duality does not necessarily extend to the thermodynamics 
of the respective universes since future event horizons exist only when a > 0. Nevertheless, 
in the particular case of pole-like expansions, Eq. ([!}, the transformation preserves the 
temperature while the entropy transforms according to Eq. (|24J) with the constraint —5/3 < 
w < — 1. More general instances will be considered elsewhere. 
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